Ferroelectric PbTiOs/SrRuOs superlattices with broken inversion symmetry 
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We have fabricated PbTiOa/SrRuOs superlattices with ultra-thin SrRuOs layers. Due to the 
superlattice geometry, the samples show a large anisotropy in their electrical resistivity, which can 
be controlled by changing the thickness of the PbTiOa layers. Therefore, along the ferroelectric 
direction, SrRuOa layers can act as dielectric, rather than metallic, elements. We show that, by 
reducing the concentration of PbTiOs, an increasingly important effect of polarization asymmetry 
due to compositional inversion symmetry breaking occurs. The results are significant as they repre- 
sent a new class of ferroelectric superlattices, with a rich and complex phase diagram. By expanding 
our set of materials we are able to introduce new behaviors that can only occur when one of the 
materials is not a perovskite titanate. Here, compositional inversion symmetry breaking in bi-color 
superlattices, due to the combined variation of A and B site ions within the superlattice, is demon- 
strated using a combination of experimental measurements and first principles density functional 
theory. 



Artificially layered perovskite oxide superlattices pro- 
vide many opportunities to develop systems with novel 
and tunable properties [H [2]. As far as ferroelectric super- 
lattices are concerned, the insulating titanium perovskite 
oxides (e.g. PbTiOa, BaTiOg, CaTiOs and SrTiOg) have 
to date been the most popular "building blocks" , but the 
need for new functionalities, particularly related to mag- 
netism, requires the use of a wider set of materials and 
a deep understanding of the new physical phenomenon 
related to interfaces. In this letter an unconventional 
approach is demonstrated: we use a material that is 
normally metallic to play the role of a dielectric, in a 
ferroelectric-dielectric superlattice. 

The much-studied compound SrRuOs provides the 
proof of concept that metallic magnetic oxides can trans- 
form into thin-film dielectric components in certain het- 
erostructures. In bulk, SrRuOa has the distorted per- 
ovskite orthorhombic Pnma structure, is metallic and 
is ferromagnetic below Tc — 160 K [3H7]. It is also a 
commonly used electrode material for oxides, and the in- 
terface with ferroelectric oxides has been much studied 
[SHU]. However, SrRuOa becomes insulating in layers of 
thickness less than 4 unit cells; this behavior has been 
observed in thin films [IMS] and in SrTiOa/SrRuOa 
superlattices [TBI ITT] . First-principles investigations [TSl - 
[20] and experiment [211 - f26] indicate that epitaxial strain, 
size effects, chemical pressure, surface reconstruction and 
interaction with the substrate may all play an impor- 
tant role in the observed behavior. In a recent con- 
tribution Vcrissimo-Alves et al. [37] showed from first 
principles calculations that a highly confined 2DEG is 
formed at the interface in SrTiOs/SrRuOa superlattices. 
We will show in this letter that a similar effect occurs in 
PbTiOs/SrRuOa superlattices, but that in the direction 
perpendicular to the interfaces PbTiOa/SrRuOs super- 
lattices containing single unit cell layers of SrRuOa are 
insulating and can be ferroelectric. 

A second motivation for creating PbTiOs/SrRuOa su- 



perlattices is that, as they have both A and B site vari- 
ation, inversion symmetry can be compositionally bro- 
ken [23 [23 • As a result, an asymmetry is introduced 
in the ferroelectric double-well potential which can lead 
to "self-poling" materials. Self-poling ferroelectric ma- 
terials are useful in piezoelectric applications where the 
desired mode of operation is to apply an electric field 
either with or against a fixed polarization direction to 
achieve, respectively, an expansion or contraction of the 
material. In a tri-color superlattice compositional break- 
ing of inversion symmetry can occur with only A or B site 
variation. The effect has been seen with A site variation 
in tri-layer superlattices containing BaTiOa , CaTiOs and 
SrTiOs [30, •' > 1 . and shown to be an appealing route to- 
wards magnetoelectric materials with tri-color variation 
on the B site [32J [321, t)ut until now this behavior has not 
been shown experimentally in bi-color superlattices. 

Epitaxial growth of PbTiOs/SrRuOs superlattices can 
be achieved on SrTiOa substrates as both PbTiOa and 
SrRuOs have in-plane lattice parameters close to that 
of SrTiOa, which at room temperature is cubic with 
a=3.905 A. At room temperature bulk PbTiOa is tetrag- 
onal (a=3.904 A, c=4.15 A) and orthorhombic SrRuOa 
can be considered as pseudo-cubic with a=3.93 A. For 
this study, the n/1 PbTiOa/SrRuOs (n unit cells (u.c.) 
PbTiOs/l u.c. SrRuOa) superlattices were deposited us- 
ing off-axis RF magnetron sputtering on (001) SrTiOa 
substrates, which had been treated with buffered HF and 
annealed to ensure Ti02 termination. The SrRuOa layer 
thickness was grown to 1 u.c. for all the samples con- 
sidered here, with the aim being that the SrRuOs layers 
should act as dielectrics, rather than metals. By con- 
trast the thickness of the PbTiOa layer was changed from 
sample to sample, so that the relative effect of bulk ferro- 
electricity vs. interfacially driven compositional inversion 
symmetry breaking could be assessed. The total number 
of bilayers in the superlattice was varied from one sample 
to another so that the total thickness of each sample was 
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in each case as close as possible to lOOnm. Growth rates 
for the two materials within the superlattice were ob- 
tained from x-ray diffraction measurements performed on 
many preliminary samples. Bottom SrRuOs electrodes 
(20 nm in thickness) were deposited in situ for the sam- 
ples used for electrical measurements and gold top elec- 
trodes were added to the samples post-deposition. The 
superlattices were grown at a temperature of 550° C and 
the SrRuOs electrodes were grown at 620° C. 

Experimentally, it has been demonstrated by Rjinders 
et al [33] that when grown by pulsed laser deposition the 
termination of a SrRuOa film is affected by both depo- 
sition conditions and layer thickness and that as a film 
grows on Ti terminated SrTiOa there is a conversion from 
a Ru02 to a SrO termination layer. At a growth temper- 
ature of 700° C this transition already occurs for a single 
unit cell SrRuOs layer, but at lower deposition temper- 
atures, this transition occurs later in the growth, as the 
Ru02 layer is comparably more stable at these condi- 
tions. As a consequence of the small thickness of our 
SrRuOa layers, the low deposition temperatures used in 
our process, and the different kinetic regime of sputtering 
as compared to pulsed laser deposition, a Ru02 termina- 
tion of our SrRuOa layers, may still be possible. 

The epitaxial growth of our samples was confirmed by 
x-ray diffraction and high resolution scanning transmis- 
sion electron microscopy (HR STEM) . Fig. [l]^b) shows a 
HR STEM cross section of an 8/1 PbTiOs/SrRuOa su- 
perlattice. The PbTiOs layers are the brightest because 
of the high atomic number of Pb, while the SrRuOs layers 
are less bright than PbTiOs, but have enhanced bright- 
ness compared to SrTiOa (not shown) because the Ru ion 
has a higher atomic number than Ti. In addition to the 
STEM image shown, we also carried out STEM-EELS 
line scans |35^ , which support the ideality of our grown 
structures. 
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FIG. 1. (a) A HR-STEM image of an 8/1 PhTiOs/SrRuOz 
superlattice. (h) The three types of interfaces considered 
theoretically in this paper, illustrated for the case of a 3/1 
PbTiOa/SrRuO^ superlattice. 

Although the interface most likely to form in 
PbTiOs/SrRuOa superlattices if the atoms in the mate- 
rial are deposited in the same ratio as the parent targets 



breaks inversion symmetry it is possible to consider in- 
terfaces that might form which do not. In Fig. [l] (b) 
we show three kinds of interfaces which we have studied 
using first principles calculations, which we illustrate for 
the case of a 3/1 superlattice. The ideal structure which 
breaks inversion symmetry is the first from the left in Fig. 
[T] (b) . In order to study the significance of the symme- 
try breaking effect, we also simulated two different kinds 
of interface which conserve symmetry, one in which one 
Pb-0 plane has been replaced by a Sr-0 plane, which 
is shown in the center of Fig. [l] (b), and another, (less 
likely due to the high volatility of Ru) , unit cell in which 
a Ti-0 plane has been replaced by a Ru-0 plane is shown 
on the right. 

We investigated the 3 kinds of interface shown above 
using first principles calculations. These were performed 
using density functional theory, using a basis of nu- 
merical atomic orbitals as implemented in the siesta 
code. We used the same basis set and pseudopoten- 
tials as Verissimo-Alves et al. [27]. We studied the in- 
fiuence of spin polarization, the use of the generalized 
gradients approximation within the commonly used Wu- 
Cohen parametrization[55] . and the effect of correlations 
within the LDA-I-U and LSDA-I-U approximations [55] . 
The different approximations used can affect the elec- 
tronic properties of the metallic layer along the parallel 
direction. However the electrical anisotropy is mostly 
dependent on the superlattice periodicity, and therefore 
results are computed using the local density approxima- 
tion (LDA), which is the best method to characterize 
both the structural and electronic properties of the su- 
perlattices. Full details regarding the calculations can 
be found in the supplemental materials to this letter [55]. 
We examined the electrical conductivity of the superlat- 
tices both in-plane ( a^x) and out-of-plane (azz) by cal- 
culating the diagonal elements of the conductivity tensor 
within the relaxation time approximation to the Boltz- 
man transport equation |40]. The anisotropy in a is fully 
determined by the anisotropy of the Fermi surface geom- 
etry, as determined by CTqq = -(^'^Tj2k'"ka^i''F-£k)- 
The relaxation time, r is the only variational parameter 
in the expression, we choose t — 1.3 x 10""'^'* s after exper- 
imental results in bulk SrRuQc. [41] . This approximation 
ignores the anisotropy of the electron-phonon scattering, 
although this is known to be a much smaller effect [10]. 
The results shown in Fig. [2] (a) are for the ideal inter- 
faces, but these quantities were also calculated for the 
two other cases and are qualitatively similar, again an 
indication that the anisotropy is almost fully determined 
by the inter- SrRuOs layer distance. We find, in agree- 
ment with ref. ^7]i that the electrons in the single unit 
cell layers of SrRuOs are confined to that layer, so while 
the in-plane conductivity, axx , of the structures does not 
change dramatically as the spacing between the layers is 
varied the out-of-plane conductivity, azz, decreases ex- 
ponentially with a characteristic length of 1.3A as the 
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thickness of the PbTiOs layers is increased. 

In the two-component PbTiOa/SrRuOs superlattices 
in which the inversion symmetry is broken by the ideal 
interface structure, our calculations predict a self-poling 
behavior. In Fig. [2] (b) we show the energy of the super- 
lattice as a function of the polarization. We computed 
the total energy of the system for atomic displacements 
along the line r = rp^ + u{fp_ — fp^), interpolating lin- 
early between the two minima of the energy. These two 
minima are characterized by two different polarization 
states, Pup (higher energy minimum) and Pdown (lower 
energy minimum) . As shown in Table |lj the simulations 
for the ideal interface show that, as the PbTiOa layer 
thickness is reduced, there is an increasingly large differ- 
ence in the values of the stable up and down polariza- 
tions, until for the 5/1 superlattice, when the potential 
well has just one minimum, only the down polarization 
is stable. 
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FIG. 2. (a) Calculated conductivity both in-plane and out of 
plane plotted a function of the number of unit cells of PbTiOz 
included in each bilayer (with zero corresponding with SrRuOs, 
strained in-plane to the SrTiO-j, lattice parameter), (b) Cal- 
culated total energy per unit cell as a function of polarization. 
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TABLE I. Calculated stable polarization magnitudes (DFT 
LDA ) in the up and down directions for a selection of ideal 
superlattices. 

The preference of one polarization state over another 
for superlattices with broken compositional inversion 
symmetry was seen in all of the calculation schemes used. 
In the spin polarized calculations for samples with single 
unit cell layers of SrRuOs the spin polarization is not 
affected by the direction of the polarization. However, 
in simulations of superlattices with symmetry breaking 
interfaces that have SrRuOa layers thicker than one unit 
cell the magnetization is different for the two polariza- 



tion directions. Although it is not the focus of the 
present paper, this finding demonstrates the potential for 
the compositional breaking of inversion symmetry at the 
PbTiOs/SrRuOa interface to enable a form of coupling 
between magnetism and ferroelectricity. 

Experimental values for the switched ferroelectric po- 
larization of the samples were obtained from polarization- 
electric field hysteresis loops performed on a number of 
samples. Polarization hysteresis was observed in samples 
with a PbTiOa layer thickness of 5 u.c or greater. The 
experimentally measured polarization as a function of the 
total number of unit cells in each bilayer is shown in Fig. 
[3j along with characteristic loops measured at 3 different 
frequencies on the 7/1 shown as an inset. Successful hys- 
teresis loops confirm that the thin layers of SrRuOa in 
the material are acting as dielectric layers, and allow a 
continuous polarization in the structure. An independent 
confirmation of ferroelectricity comes from x-ray diffrac- 
tion reciprocal space maps around superlattice Bragg 
peaks, shown in the supplemental information [3S j . These 
show diffuse scattering from the in-plane periodicity of 
stripe domains with polarization oriented up and down 
with respect to the substrate, and are similar to those 
seen in PbTiOa/SrTiOs superlattices. 02 iS] • These do- 
mains features were observed in the 7/1, 9/1, and 13/1 
superlattices, but not in the 5/1 superlattice. The lack of 
domain features in the 5/1 PbTiOa/SrRuOs superlattice 
may be due to the relative instability of one polarization 
state with respect to each other, ie., while the polar- 
ization can be switched under field from one direction to 
the other, its equilibrium configuration is dominated by a 
single polarization direction. Direct comparison between 
the calculated stable polarizations in the layer and the 
experimentally measured switched polarizations are diffi- 
cult to make, as the two quantities, while related, are not 
identical. However, it can be seen that the polarization 
values and their dependence on the number of unit cells 
in each bilayer are qualitatively similar for experiment 
and the 3 theoretical cases considered in Table lU For an 
equivalent composition the Ru excess superlattices have 
higher polarizations than the Sr excess superlattices. In 
the ideal case marked differences in the polarization value 
for the up and down state become noticeable for super- 
lattices whose bilayers contain less than 7 unit cells of 
PbTiOs. When the PbTiOs layer thickness is 3 layers or 
less the samples become fairly conductive in the out of 
plane direction, both in theory and experiment. 

An indirect probe of the stable polarization is the av- 
erage tetragonality (c/a) of ferroelectric superlattices, 
which can be measured using x-ray diffraction |44H46j . 
Our experimental measurements and first principles re- 
sults from DFT LDA calculations are shown in Fig. |4] 
As with the polarization values in Fig. [3) we have plotted 
the results in terms of the total number of unit cells per 
bilayer, and included for each set of values the nominal 
composition. However, the Sr and Ru excess structures 
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FIG. 3. Switched polarization plotted as a function of the total 
number of unit cells m each bilayer. Inset: Polarization-field 
hysteresis loops measured at 3 different frequencies on a 7/1 
sample. 



deviate from these ideal structures as described earlier; 
for the former one Pb-0 plane has been replaced by a 
Sr-0 plane, and for the latter a Ti-0 plane has been re- 
placed by a Ru-0 plane. It is difficult to make a definite 
conclusion which interface is present in our experimental 
samples solely from comparing the experimental results 
with the theoretical predictions shown. However, we sug- 
gest that a comparison of the data over the whole range 
of the plot would tend to exclude the Sr excess interface, 
and although the Ru excess interface matches the exper- 
imental data relatively well, this interface is unlikely to 
occur in experiment due to the high volatility of Ru. An 
interesting point in this figure is that the average tetrago- 
nality begins to increase again as the number of unit cells 
in each bilayer is decreased. This is because, in contrast 
to PbTiOs/SrTiOs where SrTiOs grown on SrTiOs has 
a tetragonality of 1 and the c/a montonically decreases, 
as the amount of PbTiOa decreases, SrRuOs grown on 
SrTiOS has a tetragonality of 1.03. It appears that the 
c axis lattice parameter of paraelectric PbTiOa in these 
superlattices is quite low, both in experiment and theory 
and is certainly well below the value of approximately 
1.03 that is predicted by Landau theory 03]. The up- 
turn in the c/a value occurs for the Ru excess samples 
first as these are, for each sample, essentially half a unit 
cell closer to being SrRuOa than the corresponding ideal 
superlattice. 

The effect that the compositional breaking of inversion 
symmetry has on functional properties is most evident in 
the dielectric response of the samples. In Fig.[5j we show 
dielectric constant measured as a function of electric field 
(at a frequency of lOkHz). The measurements show an 
evolution from a typical butterfly loop for PbTiOs rich 
samples (13/1 and 9/1 samples) to a highly asymmet- 
ric curve with greatly enhanced peak dielectric constant 
for the 3/1 sample. An unusual characteristic where two 
peaks are seen on each voltage trace displayed in the 5/1 
sample. The composition at which the transition from 
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FIG. 4. Average tetragonality, or the ratio of the average 
out of plane lattice parameter c to the lattice parameter of 
the SrTiOs substrate a, (i.e. c/a), plotted as a function of 
the total number of unit cells in each bilayer. As well as the 
experimental results results from DFT LDA calculations for 
three different kinds of interfaces are shown. 



conventional ferroelectric behavior occurs (approx. 7/1) 
matches the composition highlighted by theory at which 
compositional inversion symmetry breaking becomes a 
dominant factor. To summarize our picture of this sys- 
tem: The 3/1 sample can be characterized as a sponta- 
neously polarized non-switchable insulator (or, in other 
words, an interfacially driven pyroelectric) . Samples be- 
tween 7/1 and 3/1 are best described as asymmetric fer- 
roelectrics. Samples with PbTiOs layers thicker than 7 
unit cells, while still containing interfaces which compo- 
sitionally break inversion symmetry are not greatly af- 
fected by them. 
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Besides the direct results presented here, our study of- 
fers a general demonstration of the possibilities unlocked 
by expanding the material set used in ferroelectric su- 
perlattices. Our findings should motivate a broader ex- 
ploration of candidate materials for the development of 
new artificially layered ferroelectrics. In particular, be- 
sides the self-poling behavior that compositional broken 
inversion symmetry produces, the use of thin metallic 
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materials as dielectric layers has intriguing potential for 
the development of highly coupled multiferroics. 
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